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Abstract 

The sensitivity of the hghtest supersymmetric particle rehc density calculation to different 
cosmological scenarios is discussed. In particular, we investigate the effects of modifica- 
tions of the expansion rate and of the entropy content in the Early Universe. These 
effects, even with no observational consequences, can still drastically modify the relic 
density constraints on the SUSY parameter space. We suggest general parametrizations 
to evaluate such effects, and derive also constraints from Big-Bang nucleosynthesis. We 
show that using the relic density in the context of supersymmetric constraints requires a 
clear statement of the underlying cosmological model assumptions to avoid misinterpre- 
tations. On the other hand, we note that combining the relic density calculation with 
the eventual future discoveries at the LHC will hopefully shed light on the Very Early 
Universe properties. 



1 Introduction 

Cosmological observations reveal that the mass content of the Universe is mostly composed of 
dark matter of unknown, but non-baryonic nature. New physics models, such as supersymmetry 
(SUSY), provide stable particle candidates for dark matter, and one can compute their present 
energy density, the relic density |T]. This relic density is often compared to the dark matter 
density deduced from cosmological observations in order to constrain new physics parameters 
(see for example [2j). The usual assumption in doing that is that the Universe is ruled by 
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the standard model of cosmology, which assumes that radiation energy density and radiation 
entropy density dominate the Universe properties in the Very Early Universe. However, before 
Big-Bang nucleosynthesis (BBN) many phenomena could have modified the physical properties 
of the Universe, such as the expansion rate, or the entropy evolution, or even the non-thermal 
production of relic particles. The calculation of the relic density is altered in these cases, and 
for example the infiuence of a quintessence-like scalar field [3], or reheating and non thermal 
production of relic particles due to the decay of an infiaton-like scalar field |1] have been 
discussed in the literature. Scenarios involving dark fiuids or extra-dimensions modifying the 
expansion rate of the Universe have been also considered in [5j. On the other hand, with the 
start-up of the LHC, we can hope for new information on the physics beyond the standard 
model, providing more hints for the determination of the nature of dark matter [6]. 

In this paper, we consider the calculation of the relic density beyond the cosmological 
standard model and propose a generalized parametrization of the modification of the entropy 
evolution in the Early Universe. EflFects of the modification of the expansion rate was studied 
in , and we derive here the necessary equations to compute the relic density in a more generic 
way. We also discuss the importance of the LHC data in the context of cosmology and claim 
that future discoveries of the LHC can lead to a better understanding of the Universe properties 
before BBN. 

To illustrate the consequences of modifications of the cosmological model on the calculation 
of the relic density we consider in the following the minimal supersymmetric extension of the 
Standard Model (MSSM) with i?-parity conservation and show the implications on the SUSY 
parameter interpretation and constraints. 



2 Relic density calculation 

The density number of supersymmetric particles is determined by the Boltzmann equation, 
which in presence of non-thermal production of relic particles takes the form: 

^ = _3Hn - {av){n' - n%) + Nd , (1) 

where n is the number density of supersymmetric particles, {av) is the thermally averaged 
annihilation cross-section, H is the Hubble expansion rate and n^q is the supersymmetric par- 
ticle equilibrium number density. The term provides a parametrization of the non-thermal 
production of SUSY particles which is in general temperature-dependent. The expansion rate 
H is determined by the Friedmann equation: 

H' = ^iPrad + PD). (2) 

Prad is the radiation energy density, which is considered to be dominant before BBN in the 
standard cosmological model. Following 0, we introduced in Eq. ([2j) to parametrize the 
expansion rate modification, pjj can be interpreted either as an additional energy density term 
modifying the expansion {e.g. quintessence), or as an effective energy density which can account 
for other phenomena affecting the expansion rate {e.g. extra-dimensions). 
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The entropy evolution can also be altered beyond the standard cosmological model, and we 
write the entropy evolution equation in presence of entropy fluctuations: 



ds 



= -3Hs + S 



D 



(3) 



where s is the total entropy density. T>d in the above equation parametrizes effective entropy 
fluctuations due to unknown properties of the Early Universe, and is temperature-dependent. 

In the following, for the sake of generality, the three new parameters N^^ Pd and T^d are 
regarded as independent. Entropy and energy alterations are considered here as eflFective eflFects, 
which can be generated by curvature, phase transitions, extra-dimensions, or other phenomena 
in the Early Universe. In a speciflc physical scenario, these parameters may be related, as 
for example in reheating models [4J. However, the large number of unanswered questions 
in the pre-BBN epoch and the complexity of particle physics models, which involves many 
diflFerent flelds, can doubt the simplicity of reheating models. In particular many open questions 
remain in inflation, leptogenesis and baryogenesis scenarios. Therefore, a complete and realistic 
description of the Early Universe would rely on several different flelds with complementary roles, 
far beyond the decaying scalar fleld description of inflation and reheating models, and the direct 
dependence between energy and entropy would in such cases be very difficult to determine at the 
time of the relic freeze-out. For this reason, we prefer to adopt a more conservative and effective 
approach in which the effective energy and entropy densities are considered as independent. 

The radiation energy and entropy densities can be written as usual: 



Prad - 3eff{T)—T'^ , 



27r^ 3 
Srad = heff{T)—T . 



(4) 



We split the total entropy density into two parts: radiation entropy density and effective dark 
entropy density, s = Srad + sd. Using Eq. (|3]) the relation between sd and E^) can then be 
derived: 

Iheff dSD 



(5) 



with 



9. 



1/2 



h. 



eff 



1 + 



\f9ln \ 3/ie// dT 

Following the standard relic density calculation method [1], y = n/s is introduced, and Eq. 
yields 



(6) 



dY 

dx 



mi 



sp 



TT 



45G 



1/2 
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where x = misp/T, rriisp is the mass of the lightest supersymmetric relic particle, and 

sd ~ _ Pd 



sd = 



Kff{TY-iT^ 



Pd = 



(7) 



(8) 
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and 

with i running over all supersymmetric particles of mass rrii and with gi degrees of freedom. 
Sd and S^) are not independent variables and are related through Eq. In the limit where 
Pd = = = = 0, standard relations are retrieved. Integrating Eq. ([71), the relic 
density can then be calculated using: 

n/t^ = = 2.755 X lO^YomijGeY , (10) 

Pc 



where the subscript refers to the present values of the parameters. Using Eqs. flllfTOl) the relic 
density in presence of a modified expansion rate, of entropy fiuctuations and of non-thermal 
production of relic particles, can be computed provided po^ Njj and Sd or T^jj are specified. 
For pd we follow the parametrization introduced in Ref. |5]: 

PD = l^ppTad{TBBN){T /TbbnY' , (H) 

where Tbbn is the BBN temperature. Different values of Up leads to different behaviors of 
the effective density. For example, = 4 corresponds to a radiation behavior, = 6 to a 
quintessence behavior, and np > 6 to the behavior of a decaying scalar field. is the ratio of 
the effective energy density to the radiation energy density at BBN time and can be negative. 
The role of is to increase the expansion rate for p^) > 0, leading to an early decoupling and 
a higher relic density, or to decrease it for p^) < 0, leading to a late decoupling and to a smaller 
relic density. Requiring that the radiation density dominates during BBN implies \np\ <C 1. 
Moreover, > imposes \pd\ < Prad for Pd < 0, strongly limiting the interest of negative Kp 
as mentioned in 

To model the entropy perturbations, we choose to parametrize in a similar wajH: 

Sd = l^sSrad{TBBN){T /TbbnY' - (12) 

This parametrization finds its roots in the first law of thermodynamics, where energy and 
entropy are directly related and therefore the entropy parametrization can be similar to the 
energy parametrization. As for the energy density, different values of lead to different 
behaviors of the entropy density: = 3 corresponds to a radiation behavior, = 4 appears in 
dark energy models, ~ 1 in reheating models, and other values can be generated by curvature, 
scalar fields or extra-dimension effects. /^^ is the ratio of the effective entropy density to the 
radiation entropy density at BBN time and can be negative. The role of will be to increase 
the temperature at which the radiation dominates for > 0, leading to a decreased relic 
density, or to decrease this temperature for < 0, leading to an increased relic density. For 
naturalness reason, we impose that the radiation entropy density dominates at BBN time, i.e. 
\k.s\ <^ 1- Constraints on the cosmological entropy in reheating models are derived in [7J, and 
we extend here the analyses to the general parametrization f fT2l ) using BBN data. 



■'-An alternative parametrization, more similar to the one used in reheating models [4 , can be done by 
choosing T^d instead of sd- 
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Figure 1: Constraints from Yp (left) and ^H/H (right) on the effective dark energy. The 
parameter regions excluded by BBN are located above the black lines. The colors correspond 
to different values of and ^H/H. 

A general parametrization is difficult for N^. In many reheating models a scalar field decays 
into supersymmetric particles, and the non-thermal production is therefore related to the scalar 
field density. To avoid imposing ad hoc general conditions, we choose = 0. We can however 
note that the main effect of the non-thermal production is an enhancement of the final number 
of relic particles, so that it is always possible to enhance the final relic density by assuming 
non-thermal production of relic particles. 

3 BBN constraints 

In order to make a realistic analysis of the allowed cosmological scenarios, we apply the BBN 
constraints. To compute the relevant abundances of the elements, we use a version of the BBN 
abundance calculation code NUC123 |8] updated with the NACRE ^ reaction compilation and 
modified to include the parametrization of the expansion rate and entropy content of Eqs. ([2D, 
(fTTI) and (fT2l) . We consider the rather conservative bounds of [TO] : 

0.240 <Yp< 0.258 , 1.2 x 10"^ < ^H/H < 5.3 x 10"^ , (13) 

0.57 < ^H/^H < 1.52 , ^Li/H > 0.85 x 10"^° , ^Li/^Li < 0.66 , 

for the helium abundance Yp and the primordial ^H/H^ ^H/^H^ ^Li/H and ^Li/^Li ratios. The 
most constraining observables are Yp and ^H/H^ and the constraints obtained are shown in 
Fig. [1] for (/^p, np), and in Fig. [2] for (z^^, n^). 
The BBN constraints can be summarized as: 

< 10-^-^ , < 10i-2«^-6-0 , (14) 
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Figure 2: Constraints from Yp (left) and ^H/H (right) on the effective dark entropy. The 
parameter regions excluded by BBN are located above the black lines. The colors correspond 
to different values of 1^ and ^H/H. 



(15) 



Also, for consistency with the CMB observations, we impose either rip > 4 and > 3, or 
Pd = Sd = Q ior T < Tbbn-, so that and Sd do not dominate after BBN. 



4 SUSY constraints 

We now consider the influence of the modified cosmological model on the supersymmetric 
constraints. The following computations are performed with Superlso Relic v2.7 [ll,12j. Con- 
sidering the latest WMAP data [13j with an additional 10% theoretical uncertainty on the relic 
density calculation, we derive the following favored interval at 95% C.L.: 



0.088 < VtDMh^ < 0.123 . 
The older dark matter interval is also considered: 

0.1 < VtoMh^ < 0.3 . 



(16) 



(17) 



In the following, we restrict the parameters to > 4, < /^p < 1, > 3, < /^^ < 1, 
and consider a constrained MSSM scenario. To allow more flexibility in the Higgs sector, we 
focus on the Non- Universal Higgs Mass Model (NUHM), in which the parameters consist of 
the universal (non-Higgs) scalar mass at GUT scale mo, the universal gaugino mass at GUT 
scale mi/2, the trilinear soft breaking parameter at GUT scale A^^ the ratio of the two Higgs 
vacuum expectation values tan/3, the p parameter and the CP-odd Higgs mass m^. 
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Figure 3: Influence of the presence of an effective energy density (left), an effective entropy 
(center), and both an effective energy with rip = 6 and an effective entropy with = 5 (right). 
The colors correspond to different values of Qh^. The black lines delimit the regions favored by 
WMAP. The favored zones are the lower left corners for the first two plots, and between the 
black lines for the last plot. 

We consider first the NUHM test-point (mo mi/2 1 TeV, tua = M = 500 GeV, Aq 0, 
tan/3 = 40) which gives a relic density of Qh^ ^ 0.11, favored by the WMAP constraints. 

Three different effects are presented in Fig. [31 the first plot shows the infiuence of the 
presence of an additional effective density on the computed relic density. We note that when 
hip and Up increase, the relic density increases up to a factor of 10^, as already noticed in [5]. 
The second plot illustrates the effect of an additional entropy density, in absence of additional 
energy density. Here when /^^ and Ug increase, the relic density is strongly decreased down to 
a factor of 10"^"^. The third plot is an example of both additional energy density with Up = 6 
and additional entropy density with = 5. In this case, the values of the relic density varies 
from 10""^ to 10^, and we notice a narrow strip between the WMAP lines in which the entropy 
and energy effects almost cancel, leading to a degenerate zone with Qh'^ ^ 0.11. 

We then study the effects of our parametrizations while scanning over the NUHM parameter 
space. About one million random SUSY points in the NUHM parameter plane {ji^rriA) with 

= mi/2 = 1 TeV, Aq = 0, tan/3 = 40 are generated using SOFTSUSY v2.0.18 fTTJ, and 
for each point we compute fiavor physics observables, direct limits and the relic density with 
Superlso Relic v2.7. 

In Fig. m the effects of the cosmological models on the relic density constraints are demon- 
strated. The first plot is given as a reference for the standard cosmological model, showing the 
tiny strips corresponding to the regions favored by the relic density constraint. In the second 
plot, generated in a Universe with an additional energy density with K.p = 10~^ and Up = 6, the 
relic density favored strips are reduced, as already shown also in [5j, because the calculated relic 
densities are decreased in comparison to the relic densities computed in the standard scenario. 
The next plots demonstrate the infiuence of an additional entropy density compatible with 
BBN constraints. The favored strips are this time enlarged and moved towards the outside of 
the plot. This effect is due to a decrease in the relic density. 
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Figure 4: Constraints on the NUHM parameter plane {ji^mA)^ from left to right and top to 
bottom, in the standard cosmological model, in presence of a tiny energy overdensity with 
Kp = 10~^ and Tip = 6, in presence of a tiny entropy overdensity with /^^ = 10~^ and Ug = 4, 
with Kg 10~^ and rig 4, with i^g 10~^ and Ug 5, and with K.g 10""^ and Ug 5. The 
red points are excluded by the isospin asymmetry of B ^ ^*7, the gray area is excluded by 
direct collider limits, the yellow zone involves tachyonic particles, and the dark and light blue 
strips are favored by the WMAP constraints and by the older interval (fTTI) respectively. 
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Figure 5: Constraints on the NUHM parameter plane (fi^rriA)^ in presence of a tiny energy 
overdensity with tvp = 10~^^ and Up = 8 together with a tiny entropy overdensity with tvg = 10""^ 
and 77.5 = 4 on the left, and an energy overdensity with /^p = 10~^ and = 6 as well as an 
entropy overdensity with /^^ = 10~^ and = 5 on the right. The colors are as in Fig. [H 

In Fig.[5l we consider two cosmological scenarios in which both energy and entropy densities 
are present. The energy and entropy densities have opposite effects and can compensate, and 
the similarity of the plots reveals the degeneracy between the two cosmological scenarios from 
the point of view of particle physics. However, using the BBN constraints, the scenario of the 
right plot can be ruled out. 

An important consequence of this example is that if we discover that the particle physics 
scenario best in agreement with the LHC data (or future colliders) leads to a relic density in 
disagreement with the cosmological data constraints, useful information on the cosmological 
scenario may be deduced: first, it would imply that the cosmological standard model does not 
describe satisfyingly the pre-BBN Universe. Second, combining all cosmological data, and in 
particular those from BBN, it would be possible to determine physical properties of the Early 
Universe and constrain Early Universe scenarios. As such, valuable constraints on cosmological 
models can be obtained from particle colliders. 

It is important to point out that all the scenarios previously described - apart in the second 
plot of Fig. [5]- are equivalent from the point of view of the cosmological observations: there is 
no way to distinguish between them with the current cosmological data. 



5 Conclusion 

In this paper, we have argued that the calculation of the relic density is very sensitive to the 
Very Early Universe properties such as the energy or the entropy contents. For this reason, 
the interpretation of the relic density in the context of particle physics should be done very 
cautiously, as the favored parameter zones can be completely shifted in any direction by un- 
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observable cosmological phenomena during the pre-BBN era. In we aheady showed that 
the addition of energy density can falsify the use of the lower WMAP limit. In this paper 
we demonstrated that the possible presence of additional entropy also strongly questions the 
use of the upper WMAP bound. Therefore we can conclude that using the WMAP bounds 
to constrain SUSY is very model-dependent since the standard model of cosmology remains 
very uncertain and highly questionable before BBN. On the other hand, we also notice here the 
importance of the discovery of new physics beyond the particle physics Standard Model: if the 
LHC data point to a new particle physics model providing a candidate for dark matter, com- 
bining relic density calculations with cosmological data will give constraints on the pre-BBN 
era and give valuable hints on the physics of the Very Early Universe. To conclude, the relic 
density, even if not predictive for particle physics at the moment, will hopefully soon appear as 
a new tool to design more precise cosmological models and to analyze the nature of the dark 
components of the Universe. 
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